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ABSTRACT 

A control-system design method, Quadratic Optimal Cooperative Control Synthesis 
(CCS), is applied to the design of a Stability and Control Augmentation System (SC AS). 
The CCS design method is different from other design methods in that it does not require 
detailed a priori design criteria, but instead relies on an explicit optimal pilot-model to 
create desired performance. The design method, which was developed previously for fixed- 
wing aircraft, is simplified and modified for application to a Boeing Vertol CII-47 helicopter. 
Two SCAS designs are developed using the CCS design methodology. The resulting CCS 
designs are then compared with designs obtained using classical/frequency-domain methods 
and Linear Quadratic Regulator (LQR) theory in a piloted fixed-base simulation. Results 
indicate that the CCS method, with slight modifications, can be used to produce controller 
designs which compare favorably with the frequency-domain approach. 


NOMENCLATURE 

Matrix Definitions 


A 

A 

B 

B 

B„ 

B P 

C P 

C u 

C x 

D 

Fr 

F 2 

f 8 

F P 

F u 

F* 

G 

G a 

G/j 

G ff 

G m odt / ied 
G P 

G u 

G uu 

G* 

G xx 

R 


helicopter dynamics matrix 
augmented helicopter dynamics matrix 
helicopter control effectiveness matrix 
augmented helicopter control effectiveness matrix 
pilot-model helicopter control effectiveness matrix 
SCAS helicopter control effectiveness matrix 
pilot-model state observation matrix 
SCAS control measurement matrix 
SCAS state measurement matrix 
helicopter disturbance input matrix 

pilot-model control input weighting matrix (SCAS objective function) 
pilot-model control rate input weighting matrix (SCAS objective function) 
SCAS control input weighting matrix (SCAS objective function) 
pilot-model total state estimator gain matrix 
pilot-model state estimator gain matrix (on pilot control) 
pilot-model state estimator gain matrix (on aircraft state) 

LQR feedback control gain matrix 

SCAS feedback and feedforward control gain matrix 

pseudoinverse additional feedback control gain matrix 

pseudoinverse additional feedforward control gain matrix 

pseudoinverse total control gain matrix 

pilot-model feedback and feedforward control gain matrix 

pilot-model feedforward control gain matrix 

SCAS feedforward control gain matrix 

pilot-model feedback control gain matrix 

SCAS feedback control gain matrix 

control displacement weighting matrix (LQR objective function) 
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Rj control displacement weighting matrix (pilot-model objective function) 

control rate weighting matrix (pilot-model objective function) 

Q aircraft state weighting matrix (LQR objective function) 

Q„ aircraft state weighting matrix (SCAS objective function) 

Q p aircraft state weighting matrix (pilot-model objective function) 

V m pilot-model motor noise intensity matrix 

V p pilot-model observation noise intensity matrix 

W aircraft disturbance intensity matrix 

[ • ] T transpose of [ • ] 

( • ] -1 inverse of [ • ] 

[ • ]* pseudoinverse of [ • ] 

Vector Definitions 

u helicopter control input vector 

(u r = [* e ,« 0 ,Sc,M) 

u 0 SCAS control input vector 

Up pilot-model control input vector 

Up pilot-model state estimate of control vector 

v m pilot-model motor noise vector 

Vp pilot-model observation noise vector 

w aircraft input noise vector 

x helicopter state vector 

(x T = [ u, w, q, 9, v, p, <f>, r ]) 
x pilot-model estimate of helicopter state vector 

y„ SCAS measurement vector 

y p pilot-model observation vector 

Scalar Definitions 
E[ • ] expected value of [ • ] 

J LQR objective function 

J„ SCAS objective function 

J p pilot-model objective function 

p aircraft roll-rate, deg/sec 

q aircraft pitch-rate, deg/sec 

r aircraft yaw-rate, deg/sec 

u aircraft longitudinal velocity, ft/sec 

v aircraft lateral velocity, ft/sec 

w aircraft vertical velocity, ft/sec 

Greek Symbols 

(} aircraft sideslip angle, degrees 

5(t) Dirac delta function 

6 e longitudinal cyclic position, in. 

6 a lateral cyclic position, in. 

6 C collective position, in. 

6 r pedal position, in. 

A real eigenvalue 

<f> aircraft bank angle, degrees 

r n pilot-model neuromuscular time constant, sec 

9 aircraft pitch attitude, degrees 

u)„ undamped natural frequency 

£ damping ratio 
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1. INTRODUCTION 


In the past, the predominant method for Stability and Control Augmentation System 
(SCAS) design has been a frequency-domain based approach, in which suitable input-output 
characteristics are obtained by matching desired transfer functions using root-locus tech- 
niques. Although this approach has proven to be a very reliable and successful technique, 
there are disadvantages which become apparent when applied to more complex aircraft with 
higher degrees of coupling. Recently, there has also been increased interest in a more in- 
tegrated approach to control-system design, in which several automatic controllers in an 
aircraft are designed to work together instead of autonomously. Frequency-domain tech- 
niques are difficult to apply in this situation because of their single-input, single-output 
nature. 

These trends have led towards a more algorithmic, or automatic, approach to control- 
system design. One area that has received a great deal of attention in recent years has 
been Linear Quadratic Regulator (LQR) theory. Although LQR theory is well suited to the 
design of autopilots in aircraft or structural control design, it is not always adequate for 
SCAS design since the theory automatically leads to a stable closed-loop system. Successful 
SCAS design requires adequate closed-loop stability and proper pilot control response; basic 
LQR theory does not provide a direct means to meet the latter requirement. 

Another control design technique which utilizes LQR theory is model-based compensa- 
tion (ref. 1), in which LQR theory is used to minimize differences between aircraft states 
and states of a math model which possesses the desired characteristics. Although this still 
promises to be a very useful and powerful design method, there are some serious disadvan- 
tages. First, the design method requires detailed a priori knowledge of the desired aircraft 
response. Second, the structure of the control design requires on-line calculation of the 
model states for the control law implementation. Finally, like most LQR-based designs, it 
usually requires state feedback, making an estimator necessary. The combination of the 
model and the estimator leads in most cases to a very high-order control-system, requiring 
sophisticated computer equipment to do the on-line computation necessary for the controller 
implementation. 

Quadratic Optimal Cooperative Control Synthesis (CCS) is a control design method 
which also uses LQR theory, but it offers two distinct features. First, it leads to a very 
simple measurement feedback controller design. Second, it requires no detailed a priori 
design criteria, because an assumed analytic pilot-model structure is an inherent feature of 
the approach. This can be a great advantage in cases where no existing design criteria exist 
because of either the nature of aircraft being controlled, or of the task being performed. The 
need for explicit a priori design criteria is eliminated through the use of an explicit optimal 
pilot-model; all that is necessary is a description of the desired task. 

Previously, the CCS method was used to design a SCAS for the longitudinal dynamics 
of the AFTI/F-16 aircraft (ref. 2). When evaluated in a fixed-base simulation, this CCS 
design was found to compare favorably with the augmentation currently being used on the 
F-16. The purpose of the work described in this paper was to investigate the feasibility of 
using CCS methodology for control augmentation design in helicopters. To do this, it was 
necessary both to simplify and modify the design process for application to a different class 
of flight vehicle. Two different SCAS controls were designed using the CCS method on a 
helicopter; two other controllers were designed using LQR and frequency-domain methods. 
The resulting control-systems were then compared in a piloted fixed-base simulation. 

The remainder of this paper describes the CCS methodology, the experiment design and 
conduct, and the results and conclusions of the experiment. 
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2. CCS DESIGN METHOD 

For optimal control design methods, a controller is sought which will minimize some 
type of objective function. As an example, in LQR theory, a quadratic sum of plant states 
and control inputs is minimized. Similarly, in model-based compensation, a quadratic sum of 
errors between the model and plant states, and the control inputs is minimized. Cooperative 
Control Synthesis (CCS) is also based on minimizing an objective function, but it differs 
from these other two methods in what is defined to be minimized. 

The basic CCS design structure is shown in figure 1. This method incorporates an 
optimal pilot-model with the aircraft model; then a quadratic sum of aircraft states, control 
inputs, and pilot states is minimized. Another advantage of CCS over other optimal control 
design methods is the use of measurement (or output) feedback rather than full state feedback 
(for the control design). 

The CCS design method is an iterative 
procedure in which the simultaneous solution 
of both the optimal pilot-model equations and 
output feedback equations is obtained. The INPUTS 
solution method is iterative because of the use 
of output rather than state feedback; this will 
be discussed in more detail later in the paper. 

The basic design steps are described below: Figure 1 - CCS Design Structure 

1) The optimal pilot- model solution is obtained for control of the augmented aircraft (on 
the first iteration of the design, the aircraft augmentation is set to zero) 

2) The LQR output feedback solution is obtained for the total aircraft/pilot-model system. 

3) Steps 1 and 2 are repeated until a stationary solution is obtained. 

The aircraft model, optimal pilot-model, SCAS solution, and design method simplifica- 
tions and modifications are described in the following sections. 

2.1 Aircraft Model 

The linearized aircraft dynamics can be expressed in terms of the relation: 

x = Ax + BpUp + B 0 u a + Dw 

where x = aircraft states 

Up = pilot control inputs 

u 0 = aircraft controller inputs 

w = aircraft disturbance (a zero mean Gaussian white noise 

process with intensity W; E{w(t)w r (r)} = W 6[t - r)) 

and A, B p , B 0 , and D are constant matrices of appropriate dimensions. All of the vectors 
described above are perturbations about the trimmed aircraft states and control positions. 
The SCAS control law u„ and feedback measurements y„ are expressed as: 

u 0 — G„y 0 (2) 

y 0 = C*x + C„u p (3) 

where C x and C u are constant gain matrices representing the combination of aircraft states 
and control positions available for measurement and feedback, and G tt is the constant feed- 
back and feedforward gain matrix. Substituting equations (2) and (3) into equation (1) 
yields the augmented aircraft dynamics: 

x = (A + B„G a C x )x + (Bp + B a G a C„)up + Dw (4) 
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Simplifying equation (4) results in: 


x = Ax + Bu p + Dw 

where A = A + BaGaC* 

B — Bp + B a G„C u 


( 5 ) 


2.2 Optimal Pilot-Model 


Pilot modelling has been a topic of research for many years. Some of the earli- 


est work in pilot modelling consisted of modelling the pilot as a linear servo-mechanism. 


In this type of model struc- 
ture, the pilot’s behavior is at- 
tributed to that of a linear con- 
troller, with a remnant term in- 
cluded to express the portion 
of pilot behavior which is not 
consistant with the linear op- 
eration on the inputs (ref. 3). 
Some of the work based on 
this model premise has con- 
sisted of describing the pilot’s 
behavior in terms of an op- 
timal controller and state es- 
timator (refs. 4, 5, and 6). 
As described in these refer- 


MOTOR 

NOISE DISTURBANCE 



OBSERVATION 


NOISE 

<v r ) 

Figure 2 - Optimal Pilot-Model 


ences, one form of pilot-model (fig. 2) that correlates well with experimental data 
is: 


Up = G,x + G u u p + v m (6) 

x = Ax-f Bup + F a C p (x-x) + F a! Vp (7) 

Up = G*x + G u Up + F u C p (x - x) + F u v p (8) 

yp = CpX + v p (9) 

where u p = pilot’s control inputs 
y p = pilot’s observations 
x = pilot’s estimate of aircraft states 
Up = pilot’s estimate of control positions 


v p = pilot’s observation noise (a zero mean Gaussian white noise 
process of intensity Y p ; E{v p (t)v p r (r)} = V p 6(t - r)) 
v m = pilot’s motor noise (a zero mean Gaussian white noise 

process with intensity V m ; E{v m (t)v m r (r)} = V m 6(t - r)) 

C p is a constant matrix defining the pilot’s available observations, and G x , G u , F x , and F u 
are constant matrices describing the pilot’s feedback compensation and loop closures. 

Consistent with the optimal pilot-model, the pilot is assumed to control the aircraft in 
a manner which minimizes some objective function J p , expressible as: 

T 

3p ==E {riL?’ / I xT Qp x + u p Tr i u p + ujR a Up)<ftJ (10) 
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From the necessary conditions for optimality for the regulator and filter, G x , G u , F x , 
and F„ must satisfy the following equations: 


G p = [ G x G u ] = -(R J )- 1 B 1 r K 
0 = KA t + A/K - KBi(Rj) -1 Bi T K + Qx 

= PC 1 r (Vp )" 1 


F p = 


F* 

F u 


0 = PA 2 r + A 2 P - PCi r (Vp) _1 CiP + DiWiDi T 
and where Ai, A 2 , Bi, Ci, Di, Qi, and W 2 are defined as: 


( 11 ) 

( 12 ) 

(13) 

(14) 


Ai 


A B 

° 0 




Bi = [0 I] C x = [C p 0] 


Dx 


D 0 

0 I 






This pilot-model has been shown to correlate well with experimental results when Q p is 
chosen to be representative of the pilot’s task, and Rj, R 2 , Vp, and V m are adjusted to 
satisfy two conditions. The first condition is related to the human neuromuscular time 
constant and the matrix G u . Considering initially a single-input, single-output case, the 
pilot control is of the form 


« P = 9xX + 0 u «p + v m (15) 

where g u could be considered to be the negative inverse of the human neuromuscular lag 
time constant, or: 


r n = 


9» 


(16) 


Experimental correlation (ref. 4) has put r„ in the range of .10 to .25 sec. Extending 
this concept to the multi-input, multi-output case, it can be seen that the diagonal elements 
of the matrix G u should be within the range of -4 to -10 to reflect the desired motor time 
constants. R 2 can be adjusted relative to Q p to produce this result. 

The second condition concerns the pilot remnant (ref. 3), or the relationship between 
pilot control output and pilot motor noise covariances, and pilot input and observation noise 
covariances. The work of Kleinman (ref. 5) shows that the following relationships lead to an 
appropriate representation of pilot remnant: 


Vp = n Pu E{y p y p T } (17) 

V m = **E{G x x* r G/} (18) 

where Pi =.01 and p u =.003. Values for V p and V m are chosen relative to W to satisfy 
these relationships. 
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2.3 Augmentation Control Solution (SCAS) 


As was previously mentioned, the SCAS is assumed to have the form given in equations 
(2) and (3). When the pilot-model described in equations (6) through (9) is included, the 
total closed-loop aircraft-pilot system can be written in state space notation as: 


( X 

r 

1 Up 

1 A 

ii 

1 * 


Up J 

L 


A Bp 

0 

0 

i x 1 

0 G u 

Gx 

0 

J U P , 

FxCp 0 

A + B 0 G 0 C X F*C P 

Bp -|- B a G 0 C„ 

* 

Lf u C p 0 

G t - F„Cp 

G u 

Up J 


+ 


r B «i 


D 

0 

0 ■ 

0 

U. + 

0 

I 

0 

0 

0 

0 

F* 

. 0 . 


.0 

0 

fJ 

y 0 = 

[C* ' 

c u 

0 

oil 


w 

V m 


'p 


X 1 
Ur 


l Up ) 


(19) 


( 20 ) 


The feedback gain matrix G a is chosen to minimize a controller objective function J a , 
defined as: 

T 

Jo =£1^11^^ J [x r Q a x + Up r F 1 Up + ujFaUp + u 0 r F 8 u 0 ]dtJ (21) 

o 

By defining q T = [x r u p r x r uj], we can rewrite equations (19), (20), and (21) as: 


q = Aq + Bu„ + Dw 
Yo = Cq 


( 22 ) 

(23) 


(24) 


Jo = Ej r lirn^ j [q r Qq + Ua r F 3 U a ]dtj 
o 

where A, B, D, and C and w are corresponding matrices in equations (19) and (20), and: 

Q 3 

The minimizing solution u 0 is (ref. 7): 
where G a satisfies the equations: 


Q„ 0 0 01 

0 Fi + G u r FjG u G u t F,G« 0 
0 Gj^FjGu G* t FjG* 0 

L 0 0 0 0J 


Uo = G 0 y a 


(25) 


0 = 2B T KLC T + 2F s G 0 CLC r 
0 = L(A + BG 0 C) T + (A + BG„C)L + D W D T 


(26) 

( 27 ) 


where W = 


wo o' 
0 v m 0 
0 0 v p 


Since no closed form solution exists for these equations, they must be solved in an iterative 
fashion. 

The necessity of using output feedback rather than state feedback in the SCAS control 
solution can be seen by studying the equations for the optimal pilot-model ( eqs. 6 - 8 ) and 
the SCAS (eqs. 22, 23, 25). The optimal pilot-model is basically a standard LQG (Linear 
Quadratic Gaussian) regulator and estimator, in which estimation of all the aircraft states 
is required for the model implementation. If the SCAS control was also based upon state 
feedback, it would require state estimation of both the aircraft states and the optimal pilot- 
model states. On the next design iteration, the optimal pilot-model would then require 
estimation of both the aircraft states and the controller states, the dimension of which has 
been increased due to the necessity of state estimation in the SCAS. This would produce an 
increase in the number of states on every design iteration. The problem is avoided through 
the use of output feedback in the SCAS design; no additional states are introduced by the 
SCAS, and the problem remains of fixed order. 


2.4 CCS Design Variable Simplification 

The design variables in the CCS methodology are Q p , Ri, R2, W, V p , V m , Q 0 , Fi, 
F2, and Fs. Even when all of the matrices are considered to be symmetric, the application of 
CCS to an aircraft model with six degrees of freedom and four control effectors will yield 178 
individual matrix elements to be specified. At first glance, this appears to be an unreasonably 
large number of variables to use the design method on any large-order problem. However, 
constraints imposed through use of the pilot-model and other logical choices in the design 
variables lead to a far more simplified choice of design variables. This section will describe 
the simplifications which were made. 

The most direct source of simplification comes from matching parameters in the optimal 
pilot-model. As was previously noted, certain constraints must be met to apply the optimal 
pilot-model properly. The first constraint is related to the pilot neuromuscular time lag, 
r„, which has been experimentally determined to be within the range of .10 to .25 seconds. 
This can be achieved by adjusting the magnitude of the R2 matrix relative to Q p . The 
second constraint is related to the remnants of the pilot observations and control inputs; 
these constraints can be satisfied through selection of V p and V m relative to W. The last 
constraint from the pilot- model concerns the selection of the pilot state weighting matrix, 
Q p . The optimal pilot-model is only valid when Q p is chosen to reflect the task which the 
pilot is performing. The selection of Q p will be discussed later; it should now be noted that 
the choice of design variables Q p , Ri, R2, V p , and V m are effectively replaced through the 
choice of the neuromuscular time constant r„, Rj, and task description (or choice of Q p ). 

Selection of the matrices Q 0 , Ri, W, Fi, F2, and F3 is also necessary. If we assume 
the main function of a SCAS is to reduce pilot workload with limited SCAS control activity, 
a logical choice of Q a , Fi, and F 2 is: 


Qa = Qp 

(29) 

Ft = R! 

(30) 

Fj — Rj 

(31) 


This choice of matrices leads to the following objective functions: 
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( 32 ) 


J p = E ^ ^ J l xT Qp x + Up r R!U p + ujRjii p ]dt 

0 

T 

J« = J P + £1^111^^ J lu 0 T F3U 0 ]dt J (33) 

0 

As can be seen, with these choices of matrices, the aircraft objective function becomes the 
pilot objective function with an additional weighting term on the SCAS control activity. 

With these design simplifications, the unspecified design variables become the matrices 
Q p , Ri, F3, and W, along with the neuromuscular time constant r n . To study the effects 
of Ri, W, and r n , a thorough parametric study was made on the eighth-order helicopter 
model chosen for this application to determine the variation in control designs caused by 
variations in these parameters. The results indicate that the effect on the designs because 
of the parameter variations was minimal; therefore, W and Rl were chosen to be identity 
matrices, and r n was kept in the vicinity of .15 sec in the subsequent design work. 

These simplifications result in two final design matrices, Q p and F3. As was previously 
mentioned, Q p is chosen in the pilot-model to reflect the task description. In this design 
situation, F3 was varied as necessary to prevent saturation of the aircraft controls over the 
anticipated flight envelope of the aircraft. If these design matrices are further specified to 
be diagonal, then the final number of design variables drops to 12 for the problem specified 
here. 

2.5 Pseudoinverse Decoupling CCS Design Modification 

Initial design work done with the previously described CCS design method indicated 
that although the method produced good on-axis responses, it also had a tendency to produce 
some undesirable off-axis responses as well. A modification to the design process was needed 
which would preserve the on-axis responses while removing the undesirable off-axis responses. 
Many methods of altering the control design structure aimed at decoupling the final design 
were tested, but none were successful. The basic problem encountered was that every attempt 
to penalize the undesirable coupling in the objective functions led to indefinite weighting 
matrices, Q p and F3. The CCS design method requires that F3 be positive definite, and 
that Q p be positive definite or semidefinite. 

Another method was developed (similar to one used in ref. 8) which consisted of defining 
the elements in the closed-loop aircraft and control matrices responsible for the undesired 
coupling present in the design. Additional feedback and feedforward gains G/* and G // 
were produced by taking the pseudoinverse of each of the desired matrices. The equations 
are as follows: 


} 


^desired = A. actual 4“ BGy$ 

(34) 

B<ie«»red = B actual 4" BG f f 

(35) 

G /b — B ^(A<f e *i rC( j — A oc t uo |) 

(36) 

G ff ^ B^(Bj c# i rc j “ Bflctual) 

(37) 


where B* is the pseudoinverse of B, defined as: 



The modified SCAS design is then found by partitioning the SCAS control matrix G 0 into 
feedback and feedforward components G xx and G uu , and combining them with the new 
gains obtained above: 



G uu ] (unmodified CCS SCAS design) 

+ G /6 G uti + G//] (modified CCS design) 


mo 


3. EXPERIMENT DESIGN 
3.1 Simulation Vehicle/Model 


The control-system design methodologies were applied to a linearized, 6 degree of 
freedom model of the Boeing Vertol CH-47 helicopter (ref. 9). The CH-47 was chosen 
because of interest in us- 
ing a variable-stability CH- 

47 helicopter at NASA for '•* 1 ^ ^ 11111111 

eventual in-flight control- 


were lor a bU knot level- flPiHfl 

flight reference condition. 

The linear model is de- T _ TT .. 

., , . . ,. T , Figure 3 - CH-47 Helicopte 

scribed in Appendix I, and 

a summary of the open-loop eigenvalues is presented in table . 
states of translational velocity (u,v,tt;) are expressed in ft/sec atl 
in degrees, and angular rates ( p,q } r ) in deg/sec. These choic 
made to provide some rough equivalence in the weightings bet 
rates, angular rates, and attitudes. Control displacements are 
inches. 

As can be seen from ta- 
ble 1, the open-loop dynamics 
of the CH-47 includes two un- 
stable modes, the most unsta- 
ble of which is real (A = .53) 
and consists of large longitudi- 
nal velocity (u) and pitch (0) 
variations. The second unsta- 
ble mode is a complex pair (u> 

= .54, f = -.16) constituting a 
dutch-roll mode. Clearly, the 
basic aircraft is in need of some type of stability augmentation and is an excellent choice for 
design evaluation. 

3.2 Control-System Designs 


Principal Characteristics 


long, velocity and pitch 


heave-pitch 


dutch-roll 


Table 1 - OPEN-LOOP 


The task chosen for this application was the execution of an Instrument Landing System 
(ILS) approach in the presence of a disturbance. The pilot task in this case can be described 
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as limiting perturbations about the trimmed state in the presence of an external disturbance, 
which for the design case was assumed to be a white noise input to the control effectors. 

The following section describes the four design methods which were employed to design a 
SCAS for the CH-47 helicopter, and the resulting designs. The first two methods described 
are the unmodified and modified CCS methods, which will be respectively referred to as 
CCS1 and CCS2. The other two methods are Linear Quadratic Regulator (LQR) theory 
and classical/frequency-domain (FD) theory. 

3.2.1 Unmodified CCS Design (CCS1) 

In addition to Q p and F3, the matrices B 0 , W, C p , C*, and C u have to be specified 
before applying the CCS design methodology. The following assumptions were made: 

1) The pilot control input and aircraft input matrices are the same (B p = B 0 = B). This 
is a reasonable assumption since there are not any additional control effectors available 
for the SCAS design. 

2) The aircraft disturbances were in the form of control input perturbations (D = B). In 
the parametric variations done during the simplification process, it was noted that the 
resulting designs were relatively insensitive to the choice of D. 

3) Full state feedback was available for pilot observation (C p = I). The availability of all 
the aircraft states for pilot observation is reasonable under instrument conditions. 

4) All of the aircraft states and the pilot control inputs were available for feedback 
(C* r = [ I 0 ]; C u r =[ 0 I ]). This assumption was made in order to allow direct com- 
parison with the LQR control (full state feedback). 

As was indicated in the description of the CCS method, the pilot state covariance weighting 
matrix Q p must be chosen to reflect the task being performed. In this case, Q p was chosen 
to describe the execution of an ILS approach. Several iterations of choice were made, and 
evaluated by the author in piloted simulations until a Q p matrix was found which seemed 
to adequately describe the task. In all of these control designs, the matrix Fj was adjusted 
to prevent saturation of the controls within the expected range of flight conditions for the 
aircraft. Following is a description of the iterations on the choice of Q p : 

Design I - The first choice of Q p was based on the philosophy that the main desire 
of the pilot in a regulation task is to stabilize the aircraft about the trim point, 
or maintain constant pitch (0) and roll ( <j > ) attitude while minimizing sideslip (/?). 
Therefore, the first choice of Q p was made with unity weighting on 9, <f>, and f) 
(equating to a weighting value of 0.328 on v covariance). 

Design II - When Design I was evaluated in the simulator, several areas contribut- 
ing to pilot workload were identified. Airspeed control was difficult, and workload 
in the lateral directional axes was high due to oscillatory responses in both roll and 
yaw. Weightings were therefore placed on the longitudinal speed u (.02), roll-rate 
p (0.5), and yaw-rate r (0.5) covariances in the second design iteration. 

Design III - Although pilot workload was greatly reduced with Design II, oscillatory 
responses were still present in the pitch and roll axes. Also, a high degree of turn 
coordination (due to weighting on v) made small heading adjustments difficult to 
achieve. Therefore, on the next design iteration, a weighting of 0.5 was placed on 
the pitch-rate q, roll-rate r weighting was increased to 2.0, and v weighting was 
decreased to .25. 
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Final CCS1 Design - The design resulting from the weighting choices discussed above was 
considered to be a satisfactory 
control-system, since no very 
high workload areas could be 
isolated. This design was used 
for the final simulator evalua- 
tion with the other control de- 
signs. A listing of control gains 
is given in Appendix II, and 
the closed-loop eigenvalues are 
shown in table 2. 

3.2.2 Modified CCS Design (Pseudoinverse Method - CCS2) 

One problem which was encountered with the CCS1 design method was the de- 
gree of coupling (longitudinal-to-lateral and longitudinal-to-vertical) present in the closed- 
loop design. Another control-system was designed using the pseudoinverse method 
described earlier in which some of the longitudinal-to-lateral and pitch-to-heave cou- 
plings were removed. The 
“desired” aircraft matrix was 
formed from the closed-loop 
matrix by removing all the 
pitch-rate/lateral coupling el- 
ements (M^, M p , M^, and 
M r ), and the pitch-rate due 
to heave (M w ) term. The 
“desired” control matrix was 
formed by removing the pitch- 
rate/collective coupling term M$ a . The resulting control design was used for evaluation 
in the piloted simulations, and will hereafter be referred to as CCS2. The resulting con- 
trol gains are presented in Appendix II, and the closed-loop eigenvalues are shown in ta- 
ble 3. 

3.2.3 Linear Quadratic Regulator Design (LQR) 

In addition to the previously described CCS designs, it was desired to compare with a 
design obtainable using standard LQR theory. The equations used in the formulation of the 
LQR design follow. 

The aircraft dynamics are expressed as: 


A 

U n 


Principal Characteristics 

-.021 

- 

- 

longitudinal velocity ' 

-.547 

- 

- 

heave 

-1.57 

- 

- 

heave-pitch 

-2.18 

- 

- 

heave-pitch 

- 

1.09 

.61 

dutch roll 

- 

1.15 

.99 

roll 


Table 3 - CCS2 


A 

Wn 


Principal Characteristics 

.005 

- 

- 

longitudinal velocity 

-3.17 

- 

- 

heave-pitch 

- 

.98 

.59 

dutch roll 

- 

.785 

.78 

heave-pitch 

- 

1.25 

.96 

roll 


Table 2 - CCS1 


x = Ax + Bu (41) 

If full state feedback is assumed, the control which minimizes the objective function J is 
found by solving the following equations: 

0 = AK r + KA - KBR -1 B t K + Q (42) 

G = -R' 1 B t K (43) 

T 

J = E | ^lim ^ J [x r Qx + u r Ru]df| 
o 

u = Gx (44) 


96- 12 




where Q and R are weighting matrices on the state and control covariances, respec- 
tively. The LQR theory was applied in a similar manner to the CCS, in that several 
iterations of weighting matri- 
ces were tested in piloted sim- 
ulation by the author. In the 
first iteration, the weighting 
matrices Q and R were cho- 
sen to be the matrices obtained 
in the CCS1 design, Q p and 
F3. This choice was unsatis- 
factory due to saturation of all 
the controls, and several adjustments were made until a satisfactory design was achieved. 
The final LQR design was tested along with the CCS1 and CCS2 designs. The con- 
trol gains are given in Appendix II, and the closed-loop eigenvalues are shown in ta- 
ble 4. 

3.2.4 Classical/Frequency-Domain Design (FD) 

A control-system was de- 
signed by the author of refer- 
ence 10, using a priori knowl- 
edge of desirable characteris- 
tics for the defined task. The 
control gains for this design 
are shown in Appendix II, and 
the closed-loop eigenvalues are 
shown in table 5. 

4. SIMULATION EXPERIMENT 

A piloted, fixed-base simulation was conducted on a NASA Ames simulator, using a full, 
nonlinear model of the CH-47 helicopter, developed by Boeing- Vertol (ref. 9). The simulation 
facility and evaluation tasks are described in the following sections. 

4.1 Simulation Facility 

The experiment was performed on the Chair 6 fixed-base simulator at NASA Ames 
Research Center. This simulator consisted of a single-seat cockpit cab equipped with con- 
ventional helicopter controls, and a typical instrument panel (fig. 4). A 600:l-scale terrain 
board and camera visual system were presented through the cab window on a color television 
monitor. The terrain board is shown in figure 5. The simulation was also set up with an ILS 
approach (with a six-degree glideslope) to one of the runways. Outer and middle markers 
and fog were also available. A more detailed description of the ILS approach is presented 
reference 11. A turbulence model was used to simulate the effects of turbulence when desired 
during the simulation, and is further described in reference 12. 

The cockpit instrumentation consisted of a horizontal situation indicator (HSI), attitude 
director indicator (ADI), instantaneous vertical speed indicator (IVSI), and instruments to 
indicate turn and slip, airspeed, altitude (both barometric and radar), heading, and torque. 
Glideslope and localizer information were presented on both the ADI and HSI, and lights 
were provided to indicate passage over the outer and middle markers. 

4.2 Evaluation Tasks 

In addition to the ILS task, the control-systems were also evaluated in two other tasks. 
The reason for doing this was to determine the overall “robustness” of the controllers to 


A 



Principal Characteristics 

-.014 

- 

- 

longitudinal velocity 

-.57 

- 

- 

heave 

- 

1.61 

.57 

dutch roll 

- 

2.6 

.87 

roll 

- 

2.0 

.92 

heave-pitch 


Table 5 - FD 


A 

Wn 


Principal Characteristics 

-.03 

- 

- 

longitudinal velocity 

-2.11 

- 

- 

heave-pitch 

- 

.99 

.75 

heave-pitch 

- 

1.79 

.56 

dutch roll 

- 

1.37 

.95 

roll 


Table 4 - LQR 
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Figure 4 - Simulator Cockpit 


Figure 5 - Terrain Board 


varying flight conditions and tasks. Since both the LQR and CCS controllers were designed 
to prevent control saturation from occuring within the flight envelope, it was a reasonable to 
evaluate the controller designs over a wide range of tasks. The evaluation tasks are described 
in more detail below. 

Instrument Landing System (ILS) Task (figure 6a) - The ILS approach was the 
task for which the control-system was expressly designed. The approach be- 
gan in a procedure turn to establish inbound localizer tracking at a 1000- 
ft altitude. After inter- 
cepting a 6-degree glides- 
lope, the pilot tracked the lo- 
calizer and glideslope down 
to a 100-ft altitude missed- 
approach point. A missed- 
approach procedure was exe- 
cuted, consisting of a climb- 

ing left turn to establish out- Fi S ure 6a " ILS Approach Task 

bound tracking on the localizer course. The maneuver ended when the pi- 
lot had intercepted the outbound localizer course. The localizer sensitivity was 
set at 1.4 degrees for full needle deflection, glideslope at 0.7 degrees full de- 
flection. The desired performance consisted of localizer and glideslope within 
one-half full deflection when intercepted. At other times (procedure turn, 
missed approach) desired performance was: altitude control ± 100 ft, head- 
ing ± 5 degrees, and standard rate turns (3 deg/sec)± 1/2 needle deflec- 
tion. Airspeed control of 60 ± 5 knots was expected throughout the maneu- 
ver. This task was performed entirely under instrument meteorological conditions 
(IMC). 

Slalom Task (figure 6b) - The purpose of the slalom task was to expose the lateral- 
directional characteristics of the aircraft. The pilot was instructed to fly between 
the markers on a runway. The desired performance criterion consisted of ground 
path control within 10 ft, altitude at 30 dt 10 ft, and airspeed 60 db 5 knots. This 
task was performed entirely under visual meteorological conditions (VMC). 


OUTER 

MARKER 
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Figure 6b - Slalom Task 


Nap-of-the-Earth (NOE) Task (figure 6c) - The purpose of the NOE task was to 
expose the longitudinal characteristics of the aircraft. The pilot was instructed to fly 
at a constant heading over a series of 50-ft berms. Desired performance consisted of 
heading control ± 10 degrees, altitude 20 ± 10 ft except over the berms (minimizing 
exposure over berms), and airspeed 60 ± 5 knots. This was also performed under 
VMC. 



Figure 6c - Nap-of-the-Earth Task 


In all of these tasks, adequate performance was defined by doubling the desired per- 
formance tolerances. Based on these desired and adequate performance measures, the pilot 
was asked to make comments about each configuration and assign a Cooper-Harper handling 
qualities rating (ref. 13). The Cooper-Harper rating scale is shown in figure 7. In addition to 
pilot ratings and comments, time histories of control positions, controller activity, and air- 
craft parameters were made on each run. After sufficient practice on each configuration, the 
pilot was asked to fly two attempts (with a pilot option for a third) before making comments 
and assigning a rating. 


ADEQUACY FOR SELECTED TASK OR 
REQUIRED OPERATION* 



AIRCRAFT 

CHARACTERISTICS 

DEMANDS ON THE PILOT IN SELECTED 
TASK OR REQUIRED OPERATION* 1 

Excellent 
Highly desirable 

Pitot compensation not a factor for 
desired performance 

Good 

Negligible deficiencies 

Pilot compensation not a factor for 
desired performance 

Fair-Some mildly 
unpleasant deficiencies 

Minimal pilot compensation required for 
desired performance 


PILOT 


Deficiencies 

warrant 

improvement 


Minor but annoying 
deficiencies 


Desired performance requires moderate 
pitot compensation 


No ft mprovement 
I mandatory 


Major deficiencies 


Control wilt be lost during some portion of 
required operation 



Moderately objectionable 
deficiencies 

Adequate performance requires 
considerable pilot compensation 


Very objectionable but 
tolerable deficiencies 

Adequate performance requires extensive 
pilot compensation 



Deficiencies 


Maior deficiencies 

Adequate performance not attainable with 
maximum tolerable pilot compensation 
Controllability not in question 

4 


Considerable pilot compensation is required 
for control 

► require 
improvement 

1 

Major deficiencies 


Major deficiencies 

Intense pilot compensation is required to 
retain control 



'Definition of required operation involves designation of flight phase and/or 
Cooper Harper Ref NASA TND5153 suhph.ises wdh accompanying conditions. 


Figure 7 - Cooper-Harper Handling Qualities Rating Scale 
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5. RESULTS 


The evaluation maneuvers were all flown by a NASA research test pilot. A summary of 
the Cooper-Harper pilot ratings is presented in table 6. 


Augmentation 

TASK 

Type 

NOE 

Slalom 

ILS w/o turb 

ILS with turb 

OPEN-LOOP 

5 

6 

8 

* 

CCS1 

3 

4 

3 

4 

CCS2 

3 

3 

2.5 

3 

LQR 

3 

3 

3 

5 

FD 

3 

3 

3 

2 


Table 6. 


The first point to note is that the open-loop aircraft clearly has poor handling qualities, and 
that some type of SCAS is indeed necessary. The second item of interest is that, overall, 
the handling quality ratings for all of the control-systems are surprisingly similar. However, 
the pilot comments and ratings expose some noticeable differences when comparing with the 
baseline FD controller. 


The FD control produced satisfactory handling quality ratings; the pilot com- 
mented that it was easy to fly, and had excellent attitude and altitude control. 
In comparison, the LQR con- 40 1 1 20 . 

, , . r ROLL RESPONSE TO 1 INCH PITCH RESPONSE TO 

trol was also satisfactory un- lateral cyclic step i inch longitudinal 

til it was evaluated on the ILS ofd _ cyclic step 


the performance was degraded 
to adequate. The pilot com- 
mented that, overall, the re- 
sponse was somewhat sluggish 
although fairly predictable. The 


40 

ROLL RESPONSE TO 1 INCH 

20 

PITCH RESPONSE TO 


LATERAL CYCLIC STEP 

15 

1 INCH LONGITUDINAL 


□ FD 

CYCLIC STEP 

-S 20 

O LQR 

3 10 






10 

l 


5 

( 



2.5 5.0 7.5 

TIME, tec 


2.5 5.0 7.5 

TIME, tec 


although fairly predictable. The Figure 8 - Time Histories of FD and LQR Controllers 
differences in the two systems can be seen by looking at some representative time re- 
sponses. Time history responses of pitch attitude to longitudinal cyclic input and bank 
angle to lateral cyclic input for both the FD and LQR controllers are shown in fig- 
ure 8. 


As can be seen, the roll attitude to lateral cyclic responses of the two controllers are 
nearly identical. However, the pitch attitude responses to longitudinal cyclic are clearly 
different. The FD control exhibited a fairly constant steady-state response to a step input, 
while the LQR control response tended to wash out over time. Although this was apparently 
not a problem in the first three evaluation tasks, the pilot did not find it to be satisfactory 
for the ILS task when subjected to an external disturbance. 

Like the LQR controller, the CCS1 controller received ratings in the range of sat- 
isfactory to adequate, depending on the task. Performance on both the NOE and 


ILS without turbulence was con- 
sidered to be satisfactory, while 
the performance on the slalom 
and ILS with turbulence was 
only adequate. The difference in 
the two controllers characteris- 
tics are again evident in the time 
histories. The plots presented 



TIME, sec TIME, sec 


in figure 9 show the response of Figure 9 - Time Histories of FD and CCS1 Controllers 


pitch attitude to both longitudinal cyclic and lateral cyclic inputs for the FD and CCS1 
designs. 
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While the on-axis longitudinal response (pitch attitude to longitudinal cyclic) is similar, 
a large off-axis (pitch attitude to lateral cyclic) response is apparent in the CCS1 design. 
This characteristic makes the slalom task difficult to perform, since the helicopter pitches 
up and down with lateral stick inputs. This was substantiated by pilot comments on poor 
airspeed and altitude control. This off-axis response would also have been excited by the 
turbulence in the ILS task, accounting for the poor rating in that task. Although other off- 
axis effects were apparent with this controller, the most predominant was the one described 
above. 

The modified CCS design (CCS2) received comparable pilot ratings to the FD de- 
sign in all of the evaluation tasks, and similar pilot comments. One comment the pi- 
lot made specifically about the 
CCS2 design was that it ap- 
peared to have better distur- 
bance rejection capabilities than 
other designs when flown in tur- 
bulence. The similarity of rat- 
ings and comments is not sur- 
prising when looking at the time 
histories of the two controllers; 
a representative sample is pre- 
sented in figure 10. What is 
interesting to note is that the 
time histories and closed-loop 
modes are so similar, consid- 
ering they were derived using two completely different methods. The most ma- 

jor difference in the modal characteristics of the two control-systems is that the 
FD has a complex pitch mode (w = 2.0, f = .92), while CCS2 has a pair of 
real roots (A = -1.6 and -2.2) with the same basic mode shape. The other 

modes are nearly identical. The similarities in both the time responses and modal 

characteristics is even more surprising when the feedback and feedforward gains of 
the controllers are compared - they are noticeably different. The great similar- 
ity in the responses and characteristics of the two controllers lends a high de- 
gree of validation to the CCS design process. The modified CCS design process 
yielded a control very similar to the one created with existing, known design meth- 
ods. 

CONCLUSIONS 

The results of this experiment established the general validity of the CCS design pro- 
cess. With slight modifications to the CCS design process, satisfactory handling qualities 
were achieved, and a design comparable to accepted frequency-domain methods and better 
than basic LQR methods was obtained in the absence of detailed a priori design criteria. 
Demonstrated advantages of the modified CCS design method include: 



2.5 50 75 

TIME, sec 


10.0 



Figure 10 - Time Histories of FD and CCS2 Controllers 


1) Elimination of detailed a priori design criteria through use of pilot modelling 

2) Simple design method and control structure 

3) Easily achieved decoupling 

4) Ease of handling high-order models 

The CCS design method can be applied to very high order models for which traditional 
design methods become difficult to use, such as models with aeroelastic modes or rotor 
states. As such, it is well suited for multiple controller designs, such as simultaneous designs 
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of SCAS and engine controllers, which could become a very practical usage since current 
trends indicate a future emphasis on total integration of aircraft control-systems (instead of 


the current method of autonomous designs). 
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APPENDIX I 
CH-47 Linear Model 

The linearized equations are expressed in state-space form as x = Ax + Bu 
All parameters are in body-fixed coordinates. 

Velocities are expressed in ft/sec, angular rates in deg/sec, and angles in degrees. 


x= < 


( u 'i 
w 
7 

e 

V 

p 

* 


> = < 


( longitudinal velocity ^ 
vertical velocity 
pitch-rate 
pitch attitude 
lateral velocity 
roll-rate 
bank angle 


u = 


( s ‘) 


f longitudinal cyclic ' 

1 6a 

L — J 

lateral cyclic 

I 6 C 

r — 1 

1 collective 

U r J 


l pedal > 


A = 


l r ) { 

yaw-rate 

) 



’ - 0.009 

0.034 

0.046 

- 0.561 

0.001 

0.000 

- 0.076 

- 0.564 

1.757 

- 0.029 

0.002 

- 0.001 

- 0.416 

0.831 

- 1.582 

0.000 

- 0.074 

0.020 

0.000 

0.000 

1.000 

0.000 

0.000 

0.000 

- 0.001 

0.004 

- 0.002 

0.000 

- 0.072 

- 0.035 

- 0.039 

0.113 

- 0.092 

0.000 

- 0.532 

- 1.058 

0.000 

0.000 

0.000 

0.000 

0.000 

1.000 

. 0.013 

- 0.009 

- 0.057 

0.000 

0.001 

- 0.087 



* 0.048 

0.000 

0.328 

0.000 ■ 



0.843 

0.000 

- 8.151 

0.000 



23.597 

0.001 

13.394 

0.006 


r » _ 

0.000 

0.000 

0.000 

0.000 


D — 

0.941 

1.085 

0.063 

- 0.044 



1.975 

29.036 

0.011 

- 3.430 



0.000 

0.000 

0.000 

0.000 



. 1.824 

2.626 

0.262 

10 . 692 . 


0.000 

0.002 

0.000 

0.000 

0.561 

0.000 

0.000 

0.000 


- 0.001 

- 0.003 

- 0.274 

0.004 

- 1.773 

- 0.121 

0.052 

- 0.073 


Trim Condition: 


U 0 


' 100.26 

ft/sec ' 

W 0 


5.18 

ft/sec 

Q o 


0.0 

deg/sec 

©0 

Vo 

> = < 

2.96 

0.0 

degrees 

ft/sec 

P b 


0.0 

deg/sec 

$0 


-.244 

degrees 

Rq J 


. 0.0 

deg/sec > 


f 6 t0 


f -1.82 

in. ^ 

6a0 

L — . 

.197 

in. 

6,0 

* — 1 

1 4.745 

in. 

l SrQ j 


l .259 

in. J 


Maximum Control Deflections: 


j 


r ±6.5 

in. ' 


> = ◄ 

I ±4.18 

in. 

$ cmax I 


1 0 - 9.12 

in. 

^ $rmax ' 


l ±3.6 

in. > 
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APPENDIX II 
Control-System Gains 

Control Law: u = G y»x + Gy yu p 

CCS1 


CCS2 


LQR 


FD 


G /fc = 


G fb = 


G/» = 


Gfb = 


-.016 

-.006 

-.077 

-.122 

.006 

-.005 

-.005 

0.0 

0.0 

.002 

.002 

-.003 

-.046 

-.068 

.006 

-.002 

-.026 

-.042 

0.0 

.001 

.002 

. 0.0 

.001 

.003 

.004 

.028 

.009 

.020 



’ .836 

-.011 

-.087 

-. 010 - 




-.005 

.866 

-.029 

.004 



Gyy = 

-.05 

.006 

.972 

.004 




.-.005 

.006 

.005 

.889 . 



- 0.016 

-.033 

-.077 

-.122 

.004 

-.002 

-.001 

0.000 

0.002 

0.002 

0.002 

-.002 

-.046 

-.069 

.006 

-.005 

-.026 

-.042 

-.001 

.001 

.002 

.-.001 

.005 

.003 

.004 

.029 

.008 

.019 



’ .836 

-.011 

-.526 

-. 010 ' 




-.006 

.866 

.042 

0.007 



Gyy = 

-.051 

.006 

.926 

.004 




.-.005 

.006 

.072 

.888 . 



' .032 

-.031 

-.072 

-.126 

.004 

0.0 

.002 

0.0 

-.002 

0.0 

0.0 

-.036 

-.051 

-.096 

-.023 

.023 

.033 

.032 

-.006 

0.0 

-.001 

.-.002 

.002 

.004 

.002 

.116 

.011 

.038 



- o.ooo 

0.000 

0.000 

0 . 000 - 




0.000 

0.000 

0.000 

0.000 



Gyy = 

0.000 

0.000 

0.000 

0.000 




. 0.000 

0.000 

0.000 

0.000. 



’ .018 

-.036 

-.087 

-.175 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.020 

-.087 

-.140 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

0.000 

. 0.000 

0.000 

0.000 

0.000 

0.121 

0.000 

0.051 



- 1.000 

0.000 

-.568 

o . ooo - 




0.000 

1.540 

0.000 

.114 



Gyy = 

0.000 

0.000 

1.000 

0.000 




. 0.000 

-.378 
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1 . 030 . 



-. 018 ' 

-.011 

.002 

-. 091 . 


.009 ' 
-.013 
.005 
-. 095 . 
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.062 
-.001 
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